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Abstract In the framework of the "ULTIMA" project, we use ultra cold superfluid 
%e bolometers for the direct detection of single particle events, aimed for a future 
use as a dark matter detector. One parameter of the pulse shape observed after 
such an event is the thermalization time constant T^. Until now it was believed that 
this parameter only depends on geometrical factors and superfluid ^He properties, 
and that it is independent of the nature of the incident particles. In this report 
we show new results which demonstrate that a difference for muon- and neutron 
events, as well as events simulated by heater pulses exist. The possibility to use 
this difference for event discrimination in a future dark matter detector will be 
discussed. 



1 INTRODUCTION 

%e is a very appealing target matter for the direct search of non-baryonic dark 
matter''^. Various attractive features special to ^He can be pointed out: the large 
neutron capture cross section provides an inherent discrimination mechanism 
against neutrons, the particle which is normally hardest to reject. The large num- 
ber of unpaired nucleons per unity of mass gives an advantage in a large number 
of WIMP models which predict a spin dependent interaction-iS.. This sensitivity 
to the so-called axial channel makes ^He complementary to the most advanced 
detectors like Edelweiss^ and CDMS^, which are mainly sensitive to the spin in- 
dependent interaction. The ^He atoms being much lighter than the expected WIMP 
masses, the recoil energy would be limited to i? < 6 keV, allowing for an additional 
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discrimination by rejection of events above this threshold, and thus an enhanced 
signal to noise ratio. The absolute purity of %e at the working temperatures of 
«130 jUK prevents the presence of false event producing radioactive impurities 
in the sensitive medium itself. The large transparency to y-rays reduces the ex- 
pected false event rate further. Because of the existence of only one thermal bath, 
the Bogoliubov quasiparticle excitations, problems of two or more badly coupled 
thermal baths do not arise ^. 

Numerical simulations show that with a large bolometer, a sensitivity to neutrali- 
nos within a large number of supersymetric models would be obtained Addi- 
tionally those simulations show that using a matrix design of many bolometric 
cells should provide a very good overall rejection to background events. Using 
a 3-cell prototype, a certain number of the requirements of a future dark matter 
detector have already been demonstrated: A 1 keV sensitivity has been achieved**., 
with further potential of improvement by replacement of the currently used NbTi 
Vibrating Wire Resonator thermometers with specially designed, microfabricated 
Silicon resonators. Coincidences between neighboring cells, mainly produced by 
cosmic muon events, were observed and support the idea of a discrimination by 
coincident measurement. The possibility of a precise energy calibration by in- 
troduction of a well measured amount of energy using a second Vibrating Wire 
Resonator was first described in Ref.^, and many details can be found in Ref. ^ . 
This calibration is special as it does not rely on the use of a radioactive reference 
source as for most other bolometric detectors. 

For a future dark matter detector, an alternative channel of discrimination would be 
highly desirable. The main idea is the parallel measurement of either scintillation 
light or ionization charges. The ratio of energy going to scintillation/ionisation is 
observed to be different between particles interacting by electronic and nuclear re- 
coil in the case of ^He and the differences found between heat measured and to- 
tal heat release expected in superfluid ^He bolometers show that this effect persists 
in our working conditions ^. An experiment testing the realization of the ionization 
channel is in preparation. In this paper we show experimental result which indi- 
cate that a simple pulse shape analysis could provide an additional discrimination 
mechanism. 



2 EVENT SHAPE 

The heart of our particle detectors are small cylindrical copper cells of Volume 
V =0.13 cm^ filled with superfluid ^He-B, similar to one, proposed by Fisher et 
a/.ii. 

These cells are immersed inside a superfluid ^He-B bath. Thermal contact be- 
tween the cells and the bath is achieved by small orifices of S = 200 /im diame- 
ter. When a particle interacts inside the cell, the corresponding energy deposition 
will heat up the cell, and the quasiparticle density will rise. Subsequently, the cell 
will thermalize with the outer bath with a time constant T/, on the order of a few 
seconds. The main idea is to demonstrate experimentally that the discrimination 
mechanism proposed^ is feasible: strongly interacting particle, in our case mainly 
cosmic muons, which cross several cells will leave simultaneously an energy in 
all the cells along their track. On the other hand, weakly interacting particles will 
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Fig. 1 Schematic view of tiie 3 cell bolometer. The cells are filled with superfluid ■'He and 
immersed into a superfluid 'He bath, which functions as the heat sink. The thermal contact is 
obtained by the small orifices. Each cell contains a VWR thermometer and can be operated 
independently. One of the cell contains a second VWR for doing calibration pulses. The '^Co 
source in another cell was used to demonstrate the 1 keV sensitivity—. 

leave an energy in only one cell. A scheme of the 3 cell bolometer is presented in 
fig-ll] 

An important part for the functioning of the bolometers is the thermometry. 
High precision thermometry is achieved by the use of Vibrating Wire Resonators 

(VWR) '2,13 -pjjg 

we use is a fine superconducting NbTi wire, bent into a 
semi circular shape, excited by the Laplace force. When driven well below the pair 
breaking velocity, the interaction of a VWR with the superfluid ^He-B is given by 
the scattering of quasiparticles on the wire surface, and a frictional force linear 
with velocity is observed. A frequency sweep over the VWR resonance thus de- 
livers a Lorentzian lineshape. The Full Width at Half Maximum (FWHM) W{T) 
depends on the quasiparticle (QP) density which changes exponentially with tem- 
perature. In the zero field limit, W{T) writes 

WiT) = asxpi-A/kBT), (1) 

where A is the superfluid gap at zero temperature and zero field, and a a prefactor 
which depends on the geometry of the VWR and on properties of the liquid. It is 
this relation which will lead us to use in the following the terms VWR width and 
temperature interchangeably. 

After a heating event caused for example by a cosmic particle, the QP density 
inside the cell will suddenly rise, and then go back to the initial QP density by 
thermalization via the hole. This relaxation is an exponential process, with the 
time constant Zj, determined by the geometry, i.e. mainly the ratio V /S and the QP 
mean group velocity v,,. We thus can write the equilibrium VWR width after an 
event at to, of amplitude A and at base temperature Wq as 

Weq(t)=Wo+Ae~'-^0(t~to), (2) 

with 0(f) the Heavyside step function. For a flat hole, i.e. a hole with diame- 
ter ^ thickness, Tj, can be calculated, considering that effects introduced by tex- 
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Fig. 2 (Color on-line) A typical event, probably caused by a cosmic muon (circles). The data 
can be fitted (dashed line) with eq.|4] Using the obtained fit parameters A and Zb, the equilibrium 
linewidth and thus the quasiparticle density can be illustrated using eq.[2](bold line). 



tures are negligible: 

T/,=4^^ (3) 
S Vg 

As Vi, shows a r'/^ dependence, Tj, is expected to only vary slowly with tempera- 
ture. It must be pointed out that our hole has a diameter of 200 jUm and a thickness 
of 400 /im, which should increase the real value considerably due to the possibility 
of backscattering in the hole on rough surfaces, and via Andreev scattering. 

The VWR being a high Q macroscopic object, it can not immediately adapt to 
a new equilibrium position after a quick temperature change, but does so with a 
delay, determined by the time constant Tw = The measured wire linewidth is 
thus 

W^eAt)=Wo+A{e~'-^ -e-'-^)0{t-to). (4) 

A typical particle event together with a fit using this formula is presented in fig. 121 
As can be seen, this formula seems to represent very well the experimentally ob- 
tained data. One important prerequisite for the above formula to work has not 
been mentioned yet: the usual hypothesis is that the internal thermalization, i.e. 
the processes which lead from the initial, very localized, ionization to a homoge- 
neous thermal quasiparticle distribution inside the cell, is very fast compared to 
the other time constants involved. The initial very high energetic quasiparticles 
thermalize very fast by creation of new pairs of quasiparticles. The final stage of 
thermalization for near thermal QP is supposed to happen by nonlinear processes 
on the cell walls, the corresponding time constant should be on the order of the 
time of flight of the quasiparticles towards the cell walls. This time constant is 
obtained by dividing the cell dimensions by the mean group velocity and amounts 
in our case to RiO.3 ms, which is three orders of magnitudes smaller than typical 
VWR response times. As a result, a sharp slope of the rising edge should be ob- 
served which is in good agreement with experimental data (fig.|2ll. The expectation 
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is thus that internal thermalization is virtually instantaneous, and that no informa- 
tion about the nature of the influence survives long enough to have an influence 
on our measurements. 

3 PEAK SHAPE ANALYSIS 

In order to study whether all of the above is confirmed by our experimental data, 
a large number of measured peaks, acquired during about 48 h in the presence 
of an AmBe neutron source has been fitted using eq. |4] For the fits, the baseline 
Wo has been determined first by fitting the signal before an event with a constant. 
Then the main fit with eq.|4]has then been done with the three free parameters A, 
Tft, tf). As for the moment we were mainly interested in obtaining clear results for 
the Tfo parameter, the analysis is limited to relatively high energetic events with 
Zi>100 keV. An improved analysis which potentially extends the analysed region 
down to 25 keV events is a work in progress. The results as a function of tem- 
perature are presented in fig. [3] In this figure, the particle events can be separated 
by two discrimination mechanisms in three groups: Firstly, it is well established 
that neutrons have a high cross section to undergo a neutron capture reaction. This 
releases a well defined energy of 764 keV, of which about 655 keV are released as 
heat. The remaining events can be separated into events leaving or not leaving at 
the same time an energy in an adjacent cell. The events showing coincidence are 
identified as muons. The last group of events which show no coincidence corre- 
spond to either cosmic muons, crossing on its trajectory none of the other cells, or 
an event caused by other particles, most likely the elastic scattering with neutrons 
(neutron recoil event). 

The first remark which can be made is that a temperature dependence does 
indeed exist. A T^^l^ law is expected due to the temperature dependence of the 
mean group velocity, but the observed data fits better with a law. Additionally, 
the observed absolute values are for example at 140 /iK more than 5 times higher 
than the calculated value of 0.64 s using the flat hole approximation (eq.O. The 
reason for this quite strong temperature dependence might be that textures in and 
near the orifice prevent the lowest energetic excitations to leave the cell, as they 
might be forced to Andreev backscatter even by small energy barriers. 

In fig.|4]a part of the data for the temperature slice 0.140 mK< T <0.143 mK 
is plotted with the energy resolved on the y-axis. The main observation from both 
graphs is the appearance of two different bands of values. All neutron capture 
events and all heater pulses show low values of T^, with the heater pulses in av- 
erage even smaller than the neutron capture events. All muon events identifiable 
by coincident pulses show high values of T^. Non coincident low energy events 
are mainly found in the upper band, but some are found in the lower band. As 
explained, muon events can, but not necessarily have to show coincidence. On the 
other hand, the neutron source was placed in the plane horizontal to the cell align- 
ment. This, and the relatively small recoil cross section means that coincidences 
for neutron recoil events should be very rare. In the absence of a neutron source, 
the lower band completely disappears, these events with low are thus clearly 
linked to the presence of the neutron source. We thus identify in fig.|4]the region 
with E <600 keV and Tfe>3.7 s as events caused by cosmic muons, and the region 
with E <600 keV and Ti,<3.6 s as neutron recoil events. 
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Fig. 3 (Color on-line) Zj, parameter for a large number of particle events and heater pulses as 
a function of temperature. Particle events releasing an energy of about 650 keV (black squares) 
are identified as neutron capture events. Events releasing less than this energy are separated in 
events showing and not showing coincidence with the neighboring cell. Two bands of values of 
can be clearly identified, with all neutron capture events and heater pulses being in the lower 
band, and all events showing coincidence in the upper band. 
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Fig. 4 (Color on-line) For a small slice of temperature (0.140 mK< T <0.143 mK), the peak 
height is plotted as a function of the T/, parameter The calibration factor as determined by heater 
pulses is a =0.44 mHz/keV. The neutron capture events are clearly identified by the narrow band 
of energies around 655 keV (0.29 Hz). For lower energies, we identify the region of high values 
of Tft as muon events, and the region of low values of t^, as neutron recoil events. It must be 
stressed out that our energy resolution using algorithms which do not try to find good values for 
Zf, is much better (in the best cases around 1 keV) than what this graph suggests. An improved 
analysis which tries to extent this graph to lower energies is a work in progress. 
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4 Discussion and application for a future dark matter detector 



It is not easy to imagine what can be the origin of this effect. It is necessary to 
find some mechanism which lives long enough to have an influence on the order 
of seconds, and which at the same time is different for different type of particles 
depositing an energy. For instance our best guess is a delayed heat release from 
metastable triplet ^He dimers. Scintillation measurements done on scintillation in 
"^He and ^He show that their radiative life constant is about 13 s. Additionally it 
has been shown that during the rapid processes after the primary ionization a large 
amount of these dimers form. The fraction of deposited energy stored in these 
metastable state depend on the particle interacting. For light particles like muons, 
electrons and 7 this fraction is about 25%, for heavy particles it is only about 
7%. As a result, the time delayed energy release is for neutrons about three times 
smaller than for light particles. First numerical simulations based on a non radia- 
tive deexcitation, and thus a partially delayed heat release look very promising. 

The current work inscribes itself in a long term project (ULTIMA) which in- 
tends to build a dark matter detector with superfluid ^He as the target matter. One 
crucial of every dark matter detector is an efficient method to distinguish between 
events caused by ordinary particles and the weakly interacting particles which are 
searched. Numerical simulations using the GEANT4 code^ already showed that 
a very high overall rejection factor can be obtained using a high number of bolo- 
metric cells arranged in a matrix. Nevertheless, a second channel of discrimination 
would be highly desirable. The possibility to measure ionization at the same time 
as the heat release is under investigation. As the ratio heat release/ionisation is 
expected to depend on the nature of the particle interaction, this has good poten- 
tial to provide a reliable alternative discrimination, but at the cost of an increased 
complexity of the experimental setup. 

The interest for the effect presented in this paper is that it does not come at an 
additional experimental cost. Nevertheless, in order to use this pulse shape anal- 
ysis for discrimination, it must be shown in future experiments that the increased 
value of T/, exists not only for muons, but in general for particles interacting by 
electronic recoil, and that the effect prevails down to the energy range of inter- 
est (1-6 keV) which can be foreseen for neutralino dark matter. Additionally it 
is highly desirable to get a good understanding for the underlying mechanism, in 
order to optimize the future bolometers. 

And finally, the studies of time delated heat relies from exited dimers well 
support the results of experiments of "Big Bang" simulation in superfluid ^, pub- 
lished in Ref. . At that experiments the time constant of cell recovering was much 
longer, about a 60 s. and the heat relies from dimers recombination counted in the 
balance of energy as a heating. The energy losses for scintillation, suggested in 
Ref. are overestimated on factor of two or more. Consequently, the conclusions 
of article has no experimental ground. The article with proper consideration of 
energy distribution for different processes for neutrons and muons events is under 
preparation. 
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